Combinatorics, wreath products, finite space groups and magnetism 
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In this lecture some mathematical tools necessary for a proper description of the Heisenberg 
antiferromagnet are presented. We would like to point out differences between ferro- and antiferro- 
magnetic cases of Heisenberg Hamiltonian for finite spin systems. The ground-state properties are 
discussed. 



I. INTRODUCTION 



The Heisenberg model of magnetism has been investigated for years. For the finite spin system consisting of N < oo 
d • spins s we obtain the following formula 

n = - SiZj - hJ2 »' (!) 
(v) * 

where the first sum is taken over all nearest-neighbor pairs (ij) and h is an external magnetic field parallel to z-axis 
fS ,' (number of the nearest-neighbor pairs will be denoted hereafter as Af) . 

•*-> , 

A. Ferromagnet (J > 0) 

The ground state of the finite Heisenberg ferromagnet has the following properties: 

• for h — the ground state is (2Ns + l)-tuplet with total spin number S — Ns and energy per spin £^ 
~JMs 2 /N; 



• For h > the above multiplet splits into singlets and the ground-state is the one of them with maximal 

r+ 



magnetization M = Ns and the energy per site £t = £q — hs; 



• the energy gap A£ + between the ground state and the first excited state is equal to h/N, so it is proportional 
to h. 



The above properties agree with the classical results (all spins are 'parallel') and with the thermodynamic limit 
• Sh (N — > oo, h — ► 0). 

B. Antiferromagnet (J < 0) 

The classical ground state of antiferromagnet is described by the so-calledNeel state configuration \ s — s s . . . s — s). 
It is evident that this state has two-fold degeneracy. If one defines (seal) the Neel state as a state with opposite 
magnetizations in (ferromagnetic) sublattices (i.e. Sa — Sb = Ns/2 and Ma = —Mb), then the degeneracy of this 
state is Ns + 1. On the -other hand Marshall and Peirels have proved that the ground state of antiferromagnet is a 
singlet with S = M = 00. Moreover, simple calculations for finite spin systems show that the ground state has the 
following properties 

• the ground state is a singlet even for h = 0; 

• the ground state is a linear combination of all states with total magnetization M = 0; 

• the energy gap A£~ decreases for increasing h > and for sufficient large h the ground state is a state with 
S = M = 1 (for very large h in the ground state S = M = Ns — the ferromagnetic ground state is obtained) . 



1 



C. Finite Lattice Method 



Considerations of finite spin systems is veijy popular and effective method, therefore it has been frequently applied 
since the pioneer work of Bonner and FisherO. Number of states, which should be considered, is (2s + 1)^, so it grows 
very quickly (for s — 1/2 and 4x4 square lattice there are 2 16 = 65536 states!). Therefore a lot of methods are used 
in ordet,to decrease a dimensionjpf the Hamiltonian eigenproblem or to simplify the solving procedure (e.g. Lanczos 
methoda, combinatorial methodsEl and group-theoretical methodO) . In the last case a translation group of considered 
lattice is taken into account, as a rule. In contrary, we investigate also a point group and a space group (of finite 
lattice), therefore a more complete state classification scheme can be obtained!]. 

1. Example: Four spins s = 1/2 with the periodic boundary conditions 

In this case the translation group is T = C4, the point group is V — {E, a} = C s and the space group — S = C± v . 
Classification scheme for 2 4 = 16 states is given in TabJ|. £ denotes the energy per spin, S, M — total spin number and 
magnetization, respectively, and 9 £ {0_i, 60, Oi, 62}, 5 € {S ,Si} and T € {A\, A 2 , Bi, B 2 , E} are the irreps of 
C4, C s and C4U, respectively. It is worth noting that two possible decompositions of the irrep E — 0_i0Oi = So ©Si 
correspond to two choices of the basis in the irreducible subspace labelled by E. In the first case the basis is complex 
and in the second — real one, respectively. In other words all states labelled by E correspond to wave vector with 
|fc| = 1 but in the first case there is additional index — signfc = ±1, and in the second — a symmetry index a = ±1 (i.e. 
a \ Ea) = a \Ea), where a G V is a reflection). The (antiferromagnetic) ground state is given as a linear combination 

^ (1+ - +-> + 1- + -+)) - ^ (1+ + — > + 1+ - -+) + |--++> + |- + +-)) . (2) 



TABLE I. Classification of states for 4 spins 1/2, J = — 1 



e 


r 




S 


M 


8 


De; 


;eneracy 


h = 


h 




e 


A! 


Ho 








-0.50 


1 




1 




e 2 


B 1 


Ho 


1 





-0.25 


1 




1 




e 2 


B 1 


Ho 


1 


±1 


-0.25(1 ± h) 


2 




1+1 




e 2 


B 2 










0.00 


1 




1 


e_ 


1 ©6i 


E 


H ©Hi 


1 





0.00 


2 




2 


6- 


1 ©9i 


E 


H ©Hi 


1 


±1 


TO. 25/i 


4 




2+2 




e 


A! 


Ho 


2 





0.25 


1 




1 




6 


A! 


Ho 


2 


±1 


0.25(1 =F h) 


2 




1+1 




6 


Ai 


Ho 


2 


±2 


0.25(1 =F 2/i) 


2 




1+1 



2 



II. METHOD 



A. Short Description 



The most important aim of our work is to determine the ground state of (finite) Heisenberg antiferromagnet and 
its properties. It has been done for spin systems with s — 1/2 and a linear chain up to 16 spins, square 4x4 lattice, 
and a 2x2x2 cube. The results and the detail description is presented elsewhere (seed) . The main 

points of used procedure are following 

1. Find number No of states with total magnetization M = (more precisely — we calculate a dimension of 
subspace Lq containing such states); 

2. Determine this states, i.e. determine the basis B in the subspace Lq; 

3. Decompose this basis into orbits of the space group (since the space group S is a subgroup of the symmetry 
group Ejv, then one can consider the action of S^r on B)\ 

4. It can be proved that when Ns is even number then the ground state is "fully" symmetric (i.e. it transforms as 
the unit irrep), therefore from each subspace spanned on a given orbit one (the unique) such state is chosen; 

5. The eigenproblem for the operator S 2 (square of total spin S) is solved for these states (the eigenvalues of this 
operator are 0, 2, 6, ... , Ns(Ns + 1) and the equation S 2 = is the most interesting); 

6. After the above presented steps the states labelled by M = 0, r = To, and 5 = are obtained and the ground 
state is a linear combination of these states — it is determined by solution of the eigenproblem H\tp) — E \tp); 

7. As a result the ground state is obtained as a linear combination of the so-called Ising configurations and its 
properties can be easily determined (e.g. spin-spin correlations, staggered magnetization, etc.). 

It should be underlined that this procedure gives only the ground state, therefore the thermodynamics properties of 
the considered system cannot be determined. These properties can be found when one solves eigenproblems for each 
total spin number S = 0, 1, 2, ... , Ns and for each irrep of the space group. 



The first three steps can be done applying combinatorial methods. The problem is: 'Find all states with a given 
magnetization, i.e. states \m\mi . . . ra^r) fulfilling the condition = M\ Since the magnetization operator 

S z = J^i s i commutes with any a € Sjv, then the action of on the set {1, 2, . . . , N} can be considered. In the 
simplest case s = 1/2 number of states with the magnetization equal M is given by the binomial coefficient 



where k = N/2 + A/Js a number of spins with a projection m = 1/2. For s > 1/2 it can be generalized by the 
polynomial coefficientu 



where n% denotes number of spins with projection i — s. Number of states is determined by the following sum 



taking over all decompositions (no^i . . . n 2s ) with the condition n,(i — s) — M. For example, for N = 3, s = 1, 
and M = one obtains 



B. Combinatorics 
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C. Finite space groups 



The fourth step of our procedure is to determine symmetry adapted basis of subspace Lq according to the symmetry 
group of the considered Hamiltonian. It means, that only these permutations a G Sjv are taking into account which 
preserve order of points ('neighborhood). These elements form (finite) space group of finite lattice. It can be shown 
that in the one-dimensional case it is the group Cnv — Cn^C s — Dm. From it follows that for a hypercubic lattice 
in d-dimensional space the space group is given as a wreath product!] 

S = Cnv I 

where elements of the symmetric group permute axes of rf-dimensional coordinate system. The above group is 
called a complete monomial group (of degree d) of the group Cnv The considerations of a linear representations of 
this group in the subspace Lq give us the appropriate symmetry adapted basis, which is used in the next steps of our 
procedure (5-7). These steps are performed using numerical methods (solution of eigenproblem for real symmetric 
matrix) . 



III. FINAL REMARKS 



Using the above described method we have obtained, e.g., energy per site for 16-spin linear chain. Its value - 
—0.446 — is very close to the exact result in the thermodynamic limit — 1/4 — In 2 w —0.4432. Therefore, one can 
say that in the one-dimensional case a system of 20 spins is quite good approximation of the infinite system. 

Application of the proposed procedure to two- and three-dimensional spin systems requires much more intensive 
investigations of group action of the symmetric group Sat on the basis states (the Ising configurations) \m\ . . . mjv). 
It should enable to consider system of 20x20 or 10 x 10 x 10 spins for any value of spin number s. 
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